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ABSTRACT 

Aims. We present new radio observations of two regions of the spiral galaxy NGC 3627, including new radio continuum observations at 1.4 
, GHz with the Very Large Array, and also new observations in the CO line, taken with the Plateau de Bure interferometer Comparing these 
' observations with archival Spitzer and data we study the correlation of the radio continuum (RC), infrared-8yum, and CO emissions at a 

spatial resolution of 100 pc. 

Methods. We compare the point-by-point variations of the RC, CO, and brightnesses in two distinct regions of 2 kpc x 2 kpc in size of 
NGC 3627. We also present a three-dimensional fit of the observed data. 

Results. We examined scale much lower than the electron diffusion scale, where a breakdown of the correlations would be expected. However, 
I no evidence for such correlation breakdown has been found. The RC emission follows the distribution of CO well, and the widths of several 

structures, measured along slices across them, are comparable. Furthermore, we found that down to a spatial scale of 100 pc, the radio continuum 
, emission is correlated with dust emissions at 4.5, 5.8, and 8 //m, which trace different dust temperatures. We present a new perspective, a three- 
■ dimensional representation, with which to visualize and study the RC-CO-24yum and RC-CO-8/jm correlations. We fit the observed data with a 

three-dimensional line, obtaining a rms of 0.25 dex. 

Conclusions. The observed correlation enhances the complexity of the electrons diffusion, losses, and injection mechanisms, and of their 
connection to star formation processes described by molecular and dust emissions. We plan to further investigate this connection using spatially 
resolved spectral index studies at low radio frequencies where the thermal emission component is seemingly negligible. 
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1. Introduction 

The star formation process powers the emissions of spiral 
galaxies across the whole electromagnetic spectrum. Star- 
forming galaxies are luminous sources of X-ray emission due 
to close accreting binaries with a compact companion, young 
supernova remnants, and hot plasma associate d to star-forming 
regions and galactic winds dFabbianol 1989t) . Young massive 
stars are strong sources of ultraviolet (UV) and, via the ion- 
ization of the interstel lar medium (ISM), of emissions 
( lRownd& Young|l999h . The dust heated by the interstellar ra- 
diation field produces the far infrared (FIR) luminosity. At mil- 
limeter wavelengths the CO emission traces the bulk of molec- 
ular g as in which star formation occurs (e.g. jYoung & Scoville 
199 ll) . the HCN emission is also a reliable tracer of dense 
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molecular regions. The radio emission, at centimeter wave- 
lengths, from galaxies consists of a mixture of non-thermal and 
thermal component, produced via synchrotron emission and 
free-free emission from HII regions, respectively. Correlations 
between global luminosities, integrated over the entire galax- 
ies, at these different wave lengths have been re ported by sev- 
eral authors. In particular, iRanalli et al. I ^003) studied a X- 
ray/radio correlation in star-forming galaxies;,Kennicutt ( 1 9981) 
established the global Schmid t law in galaxies using t he cor- 
relation between IR and CO; iRownd & Young (Il999h found 
a strong correlation between the Hq. surface brightnesses and 
the CO integrated intensity; and, finally, a tight correlation 
between IR and HCN luminosities has been observed by 
Gao&Soloinonl(l2004 . 



The tightest correlation observed in spiral galaxies is 
that between th e non-thermal radio and the FIR luminosity 
dCondon 1992 ). This correlation extends over five decades 
and h olds for a rema rkably wide variety of galaxy types 
(e.g., lYun et alj l200lh . It applies not only globally, but 
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Fig. 1. Overlay of the VLA 1 .4 GHz (contours) and BIMA CO 
(grey-scale) intensity. The synthesized beam of both observa- 
tions (6") is shown in lower left corner. The circles mark the 
north and south bar end regions we observed with the Plateau 
de Bure Interferometer and with VLA; their diameter repre- 
sents the primary beam of a single PdBI antenna at 115 GHz. 



also on kpc scales within the disc s of individual galax 



ie s (e.g., B oulan ser & Perault) 


19881 


Bicav&Helou 19901 


Murphvetal. 2006a; Tabatabaei et al. 


2007b. It has been 


established also that the CO 


emission in spiral galaxies 



is well correlated with the radio-continuum (RC) emis 



sion o n global and intermediate scales (e.g., iRickard et al.. 
19771; llsrael & Rowan-Ro binsonI [l984t lAdleretal.1 ll99U 



Murgiaetal.1 l2002h Recently, iMurgia et all (120051) and 
Paladino et aP (|2006|) extended the study of the RC-CO cor- 
relation down to angular scale of ~ 6" for a sample of 22 
spiral galaxies selected from the Berkeley-lUinois-Maryland 
Associ ation 3mm CO Survey of Nearby Galaxies (BIMA 
SONG; Rega n et al.ll200ll; iHelfer et alj|2003h . They found that 
the RC-CO correlation is as tight as the global one, down to 
sub-kpc size scales. 

Despite the fact that the aforementioned local RC and CO 
studies on galaxies enhance differences between the slopes of 
the correlation from galaxy to galaxy, and inside the same 
galaxy, from region to region, they clearly indicate that the 
RC brightness is spatially correlated to the CO emission. So 
far, the scarce angular resolution of the FIR observations pre- 
vents detailed investigations of the FIR-RC correlation on sub- 
arcminute scales. Nonetheless, recent high resolution Spitzer 
observations allow us the comparison of RC and CO emissions 
to shorter wavelengths than the FIR ones. In fact, using Spitzer 
24 jum im ages, we found that RC is still correlated with this 
emission dPaladino et al. I I2OO6I) . The question that we are try- 



ing to address is, what are the physical basis of these spatial 
correlations. The non-thermal RC is, for the most part, syn- 
chrotron emission that arises from the interaction of relativistic 
electrons with the ambient magnetic field in which they dif- 
fuse. These electrons are supposed to be accelerated in super- 
nova remnants, which represent the ultimate phase of the life 
of massive stars. The reason why RC emission is so closely 
related to IR and CO is not so straightforward. We know that 
the relativistic electrons diffuse from their birthplaces, there- 
fore one would expect that the spatial correlations break down 
below the characteristic diffusion scale-length of the radiat- 
ing electrons. A typical value of relativistic electrons diffusion 
scale in spiral galaxies is of the order of a few kpc (see also 
Sect.|4]i. Although, the value of the diffusion scale is extremely 
uncertain, depending on the propagation of cosmic-ray elec- 
trons along magnetic field lines, which is not fully understood. 

The recent result o f a tight RC-CO-24 fim correlations 
down to sub-kpc scale dPaladino et al. 2006 ^ indicates either 
that we have not yet probed the spatial scale at which the cor- 
relations break down or that there is a mechanism of regula- 
tion that compensates the electrons diffusion. With the aim of 
testing the former possibility we prese nt the results o btained 
on one galaxy of the studied sample ( Paladino et al.l I2OO6I) . 
NGC 3627, observed at resolution of ~ 100 pc, by extending 
the study of the correlations between RC and CO emission to 
higher linear resolutions. 

We present new RC and CO observations, at a resolution 
of 2", of two regions of the spiral galaxy NGC 3627, obtained 
with the Very Large Arra}0 (VLA) and the Plateau de Bure 
Interferometer (PdBI), respectively. The galaxy NGC 3627 is 
a promising target because it is a relatively nearby , face-on 
galaxy, at the distance of 11.1 Mpc dSaha et al.lll999l) . At this 
distance 1" corresponds to 53 pc. Both RC and CO emissions 
are localized in a narrow bar structure of ~300 pc in width. The 
emissions show a peak at the position of the nucleus, which 
extend along the leading edges of the bar forming two broad 
peaks at the bar end s, and then trail-off into the spiral arms. 
Paladino et all ( I2OO6I) found that 24 /im-RC-CO correlations in 
this galaxy persist down to a spatial scales of 350 pc. The new 
RC and CO observations are focused only on two regions of 
the bar far away from the nucleus of the galaxy (see Fig. [ij. 
As NGC 3627 is a Seyfert galaxy, we prefer to avoid the core 
whose RC emission is likely contaminated by the active galac- 
tic nuclei. 

We report on the outcomes of new RC and CO observations 
of two regions of NGC 3627, at resolution of ~2", correspond- 
ing to ~ 100 pc at a distance of 1 1 Mpc. The resolution achieved 
is comparable to that of mid-IR Spitzer observations. In order to 
add information about the connection between the RC emission 
and other thermal emissions on spatial scale of 100 pc, we com- 
pare our new observations to archival 8, 5.8, 4.5 fim Spitzer and 
H„ observations. This work is organized as follows: in Sect. 2, 
we summarize new and archival observations, and give details 
of data reduction procedures. Data analysis and results are pre- 
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Fig. 2. North bar end of NGC 3627. Upper panel: '^CO(l-O) velocity-channel maps observed with the PdBI. The angular 
resolution is of 2'.'72xl'.'75 at PA=31°. We map an area of view of 35"(the diameter of the primary beam at 1 15 GHz is of 43")- 
The phase tracking center is indicated by a cross at Q;y2000=l l''20"'13.5* and (J720()0=13°00'17'.'7. Velocity-channels are displayed 
from v=-47.28 km s"' to v=-180.88 km s"' in steps of 3.26 km s"'. Velocities are referred to the LSR and the zero velocity 
corresponds to vo=712.6 km s"'. Contour levels start from -30 mJy beam"' in step of 30 mJy beam"'. The rms noise is 6 mJy 
beam"' and only regions whose brightness is larger than 5-cr are shown. 

Lower panel: North bar end of NGC 3627: '^CO(l-O) integrated intensity. The left panel shows the velocity integrated image. 
Contour levels start from 0.4 Jy km s"' beam"' (the l-cr level is 0.13 Jy km s"' beam"') and scale by a factor of V2. The right 
panel shows the image obtained using the smooth-and-mask technique (see Sect |2.1t . The grey-scale image shows emissions 
larger than the pixel 3 cr level. The contours start from 0.3 Jy km s"' beam"', an average value of the rms in fainter regions, and 
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Fig. 3. South bar end of NGC 3627. Upper panel: '^CO(l-O) velocity-channel maps observed with the PdBI. The angular 
resolution of 2'.'75xl'.'73 at PA=31°. We map an area of view of 25", i.e., about half of the diameter of the primary beam at 
115 GHz. The phase tracking center is indicated by a cross at Q'y20oo=l l''20'"16.6' and 5y2000=12°58'44'.'5. Velocity-channels 
are displayed from v=216.68 km s ' to v=60.26 km s ' in steps of 3.26 km s Velocities are referred to the LSR and the zero 
velocity corresponds to vo=712.6 km s"'. Contour levels start from -30 mJy beam"' with step of 30 mJy beam"'. The rms noise 
is 7 mJy beam"' and only regions whose brightness is above 4-cr levels are shown. 

Lower panel: '^CO(l-O) integrated intensity. The left panel shows the velocity integrated image. Contour levels start from 0.5 Jy 
km s"' beam"' (the l-cr level is 0.16 Jy km s"' beam"') and scale by a factor of V2. The right panel shows the image obtained 
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sented in Sect. 3. The results are discussed in Sect. 4, and a 
summary is given in Sect. 5. 

2. Observations and Data Reduction 

2.1. CO data 

The observations of the CO emission from NGC 3627 were car- 
ried out with the IRAM interferometer on March 20, 2005, us- 
ing the 6Cp configuration of the array. We observed simulta- 
neously the J=l-0 and J=2-l lines of '^CO at the bar ends of 
the galaxy (for details see Fig. [T] and Tab. [U; the beam size 
is 2'.'7 X 1'.'7 and 1'.'7 x a.' 8 in the CO(l-O) and CO(2-l) line, 
respectively. 

The spectral correlator was split in two IF centered at 
1 14.997 and 229.99 GHz, respectively. These are the transition 
rest frequencies of the two '^CO lines for an assumed reces- 
sion velocity of vo(LSR)=712.6 km s"', which is the average 
of the recession velocities at the galaxy bar ends (Tab.[T]). Using 
four partly overlapping 160 MHz-wide units, the correlator was 
configured to cover a bandwidth of 580 MHz for each line, cor- 
responding to 1512 km s"' and 756 km s"' at 1 15 GHz and 230 
GHz, respectively. We obtained visibilities with on-source inte- 
gration times of 10 minutes framed by short (~ 2 min) integra- 
tions on the nearby quasars used as amplitude (3C84, 3C273) 
and phase (0923+392, 1055+018) calibrators. Total observing 
times are reported in table |2l 

We are primarily interested in the CO(l-O) transition, and 
combine the PdBI data with the available BIMA observations. 
The '^CO(2-l) observations are heavily affected by the "miss- 
ing flux" problem (see Sect. 12.1. Il l, which makes them not re- 
liable for a more quantitative analysis. However, for the sake 
of completeness, the CO(2-l) velocity channel maps and inte- 
grated intensity images are also reported in Appendix A. 

In Figs.|2]and[3](upper panels), we show the velocity chan- 
nel maps observed with PdBI in the north and south bar ends, 
correspondingly. The rms derived from emission-free channels 
of 3.26 km s ' is of 6 mJy beam ' in the north bar end and of 7 
mJy beam"^ in the south bar end. No continuum emission was 
detected, down to 1 cr rms noise level of 0.4 mJy beam ' (0.5 
mJy beam '), for the north (south) bar end. 

In Figs. |2] and [3] (lower panels), we present the '^CO inte- 
grated intensity images for the two observed regions. The left 
panel shows the integrated intensity obtained by summing all 
channels in the data cube. The 1-cr level is given by the noise 
estimated in the free channel multiplied by VA^Av, where N 
is the number of channels summed and the Av is the width of 
each channel (3.26 km s '): in the north and south regions, 
the 1-cr levels are 0.13 Jy beam"' km s ' and 0.16 Jy beam ' 
km s~', respectively. The right panel shows the integrated in- 
tensity forrned us ing the smooth-and-masked technique (see 
Heifer et al l2003h . which is very efi'ective in showing low level 
emission in the map. First, we smoothed the data cube by a 
Gaussian with FWHM=6" and we created a mask accepting 
all pixels in each channel map where the signal was stronger 
than 3 cr rms noise level of the smoothed channel. Second, 
we selected in the full-resolution data cube only those pixels 
that with non-zero values in the mask file and summed them 



Source RA 


DEC 


LSR Velocity 


(J2000.0) 


(J2000.0) 


(km/s) 


NGC3627N 11:20:13.5 


+ 13:00:17.7 


+600 


NGC3627S 11:20:16.6 


+ 12:58:44.5 


+825 



Table 1. Pointing positions and local standard of rest velocities 
for the two bar ends. 





VLA* 


PdBF 


cri.4 


C CO(l-O) 




hr 


hr 


yuJy bm"' 


Jy bm ' km s ' 


NGC3627N 


6 


8 


22 


0.13 


NGC3627S 


6 


8 


24 


0.16 



Table 2. Observations summary. Observing times and rms lev- 
els for VLA and PdBI observations (see Sect l2.1l and l2.2l ). 



' Observing times include overheads. 



to obtain the integrated intensity image. Finally, we multiplied 
the pixel values by the velocity width of an individual channel 
(3.26 km s"') to generate the image in the units of Jy beam"' 
km s"'. This technique implies a variation of the pixel -by-pixel 
rms in the CO integrated intensity image. The plots in Figs.|2] 
and|3]show contours starting from 0.3 Jy beam"' km s"', the 
average value of the rms in fainter regions, whereas the grey- 
scale images show emissions above 3 times the rms measured 
pixel-by-pixel. 
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Fig. 4. Comparison of BIMA and PdBI spectra measured in the 
northern emission region. Data have been averaged in RA-DEC 
plane, obtaining the profile along the velocity axis. The PdBI 
spectrum is shown in (red) dotted line, the BIMA in (black) 
short dashed line and the short-spacing corrected spectrum in 
(blue) solid Une. 



6 



R. Paladino et al.: Correlation of the radio continuum, infrared, and CO molecular emissions in NGC 3627 




112018 112014 IIEnia 11 JQ19 11201B 112017 112018 113015 

RIGHT ASCENS ON fJ200D) RIGHT ASCENS CM fJ2CDD) 



Fig. 5. VLA-A array images of north bar end {Left) and south bar end (Right) of NGC 3627. The angular resolution is of 1 '.'9 x 
1 '.'3 at PA=49°. Contours level start from 72 /iJy beam ' and 66 //Jy beam ' (3cr), respectively, and scale by a factor of V2. 



2.1 .1 . BIMA and PdBI data combination 

Our PdBI observations are affected by the "missing flux" prob- 
lem. The baselines range from 16 to 330 meters. The observed 
regions have extension of ~ 25", larger than half the primary 
beam size of PdBI antennae at 115 GHz (43"). In this case, 
the PdBI observati ons may have missed about 50% of the flux 
(iHelfer et al.ll2002l) . 

To reduce this problem we combined the PdBI data with 
those of the BIMA SONG. The BIMA interferometer has 
smaller antennas than PdBI (dish diameter is 6 m), and min- 
imum dish separation of 8 m, the largest angular size de- 
tectable at 115 GHz being ~67". In addition, data obtained 
with NRAO 12 m telescope have been incorporated to inter- 
ferometric maps in BIMA SON G images, therefo re these im- 
ages are not missing any flux ( Heifer et al.ll2003h . The com- 
bination should be performed by merging PdBI and BIMA u- 
V data together and then c ombining thi s with the 12 m data 
for the deconvolution (see iHelfer et al.l (2003) for the appli- 
cation of this method to the BIMA SONG). In our case, us- 
ing this method is very difficult, since we have two discrete 
PdBI fields across a region containing extended emission ob- 
served in a multi-pointings mode with BIMA. We decided to 
combine PdBI and BIMA data in the image plane, using the 
technique implemented in the MIRIAD (Multichannel Image 
Reconstruction Im age Analysis and Display) task IMMERGE 
dSau lt et al. 1995). This technique requires caution, neverthe- 
less, iWong_&.£lLtz (2000) show that this method agrees well 
with the "linear combination" described in Stanimirovic et al 



(Il999h m which data are combined before deconvolution and 
then deconvolved with a composite beam. The task IMMERGE 
appUes the "Fourier" technique; the maps are transformed to 



the Fourier plane and merged together. In a specified annulus 
of the Fourier domain, where both images should agree, we cal- 
culated the flux calibration factor, which minimize differences 
between the data of the two resolutions. Since the ranges of 
baseline length of the PdB and BIMA interferometer are 20- 
230 m and 8-87 m, respectively, the annulus used is between 
30 and 80 meters. Since the primary beams of the two images 
are different, the low-resolution image is tapered to match any 
residual primary beam response in the high-resolution one. 



For the combination of PdB and BIMA data, we regrid the 
PdBI cubes to the same velocity grid as the BIMA+12m cubes. 
For comparison between the PdBI and BIMA data, in Fig.|4]we 
show the average spectrum measured in the northern region. 
Data have been averaged in RA-DEC plane, obtaining the pro- 
file along the velocity axis. In this plot it is evident that the in- 
clusion of the short-spacing data has a noticeable impact on the 
flux of channels whose velocities lie between -120 and -180 
km s"', i.e., those with more diffuse emission. Missing flux still 
remains an issue in our CO combined images. Nevertheless, 
we believe that the adopted combination technique is sufficient 
for our purpose of analy zing the moment-0 images. In partic- 
ular, iHelfeLglalJ (120031 their Fig. 53) indicate that the flux is 
uniformely lost by BIMA throughout the galaxy. Hence, if a 
similar behavior can be expected for our combined maps, the 
missing flux would only produce a horizontal shift of the point- 
to-point RC-CO correlations shown below (see Figs.|7]and[8]l, 
without affecting the slopes. Such an effect would not invali- 
date the overall results. 
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Fig. 6. Comparison between Ii 4, Ico and S/^m mid-IR emissions. Brightnesses are measured by averaging all pixel values within 
circles of 18"in diameter recovering the entire galaxy. Filled blue triangles represent the measures in the south bar end, filled red 
circles are those in the north bar end. 



2.2. Radio continuum data 

We observed the 1.4 GHz emission from the two regions of 
NGC 3627 with the VLA in its A configuration on April, 1 1 , 
2006, spending ~ 3.5 hours on each target . We observed with 
a bandwidth of 50 MHz for the two IF channels centered at 
1.465 and 1.385 GHz. 

The data were calibrated and imaged with the NRAO pack- 
age AIPS following standard procedures: Fourier- Transform, 
Clean and Restore. Self-calibration was applied to minimize 
the effects of p hase variations. Th e flux densities were brought 
to the scale of Baars et aTl ( 1977 ) using the calibrators 3C138 
and 3C286. The final cleaned maps were restored with beams 
of l'.'9xl'.'3. Figure |5] shows the VLA- A array maps obtained 
for the two observed regions, produced using natural weight- 
ing. The final noise of total intensity images obtained for the 
north and south bar end region is of 22 and 24 fiiy beam 
respectively. 

To improve the sensitivity to large scale structures, the A- 
array data were combined with D and B-array data in the visi- 
bility plane, using the A IPS task DBCON. Th e B-array obser- 
vations are described in Paladino et al. ( 20061) . while those of 
the D-array are archival data. The combined image has a reso- 
lution of 2 '.'07 X 1 '.'43 and a rms level of 21 //Jy beam ' . 

For the comparison with CO and IR images, we convolved 
the A-i-B-i-D-array VLA images to the beam of both BIMA and 
PdBI observations, i.e., 7'.'3 x6"and 2'.'7 xr.'8, respectively. 
The derived low and high resolution images have a 1-cr noise 
of 136 yLfJy beam"' and 27 //Jy beam"', correspondingly. 



2.3. Arciiival data 

The spiral galaxy NGC 3627 has been observed as part 
of the Spitzer Infra red Nearby Galaxies Survey (SINGS 
Kennicutt et alj|2003h . and IR Spitzer images. In addition, H„ 
images have been released in the context of the SINGS Second 
Data DeHvery, April 2005. 



The observations obtained with the Infrared Array Camera 
(IRAC: lFazio et al.l2004) at 4.5, 5.8, and 8 yum have resolutions 



of 1 '.'2, 1 '.'5, and 1 '.'8, respectively. The Ho- images, obtained at 
Kitt Peak National Observatory (KPNO) 2.1m telescope, using 
a narrow-band filter centered at 6618 A, have a PSF of 1 '.'9. For 
details on Spitzer and Ha observations and data processing see 
delivery document for SINGS Second Data DeliverjU. 

The resolution of IRAC and Ha images is only slightly 
smaller than that of our RC and CO images, making the com- 
parison between these sets of images consistent. To this end, we 
convolved the IRAC and Ha images to the largest beam and 
performed the geometric transformation of images to match 
with the pointings of RC and CO observations. 

We compare the RC and CO emission with the dust emis- 
sion at mid-infrared wavelengths. The images at A24pm can 
be considered pure 'dust' images, with negligible contributions 
from the photospheric emission of stars and from nebular emis- 
sion, whereas the IRAC wavelenghts (3.6, 4.5, 5.8 and 8 jum) 
are dominated by stellar emission. In particular, the 4.5 fim 
emission is dominated by photospheric emission from stars, 
but can also contain a component of hot dust emission. In 5.8 
and 8 yum images, however, the emission from PAH becomes 
significant. In order to obtain the distribution of dust at these 
wavelen gths by subtracting t he stellar contribution, we used the 
recipe of iHelou et al. i iEool . These authors assume that all the 
3.6 pm emission is due to stars and extrapolate this component 
to longer wavelengths using stellar population mod eling. Tore 



produce the star contribution at 4.5, 5.8, and 8 fum, iHelou et al 



(120041) scale the 3.6 pm map by 0.596, 0.399, and 0.232, re- 
spectively, and subtract these scaled maps from the observed 
maps pixel by pixel to yield the "dust maps". Since this method 
may underestimate the dust component e.g., the 3.6 pm images 
can contain a hot dust component in addition to photospheric 
emission we judge this approach as a reliable one for our pur- 



- http://data.spitzer.caltech.edu/popular/sings/20050506_enhanced 
_vl /Doc uments/sings_second_delivery _v2.pdf 
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poses. In fact, the hot dust emis sion has an irnpact b elow few 
percent on the dust-only images JCalzetti et alJl2005 ). 
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3. Data analysis and results 

We extend the study of the RC-CO-24 //m correlations obtai ned 
for a spatial scale of 6", reported in Paladino et al. ( 2006f) . to 
the 8 jum Spitzer data, convolving 8 fim images to the BIMA 
resolution: 7 '.'3 x 6". To obtain a general view of the correla- 
tions across the entire galaxy, we overlaid a grid of circles with 
a diameter of 18" across the spiral structure. Then, we calcu- 
lated the 1.4 GHz RG brightness (I1.4), the GO integrated in- 
tensity (Ico) and the mid-IR brightness (Ig^m) by averaging all 
pixel values within each circle. The correlations obtained are 
shown in Fig.|6] We fitted the observed data with a power law 
parametrized as y(x) - c- jc*'"'"", taking into account the errors of 
both X and y. The results of this fit and the Pearson coefficient 
of the correlations are reported as labels on the plots shown 
in Fig.|6l The correlations of the 8//m emission with both RG 



and CO have scatters 0.08 and 0.15 dex, respectively. This val- 
ues are even lower than that of RC-GO correlation. In partic- 
ular the RG-8 jum correlation has a Pearson coefficient near 1 . 
The slopes of the correlations are very different, depending on 
which variable is associated to the x and y-axis. To provide a 
more complete description of the RG-GO-8 /urn correlation, we 
introduce a three-dimensional view of it in Sect. 13. 31 

In Figs. |6] points corresponding to northern and southern 
regions, observed at high resolution with the PdBI and the VLA 
in this work, are marked with filled blue triangles and filled red 
circles, respectively. In the following section, we report on the 
correlations studied at 2" resolution in these two regions. 

3.1. The mid IR-RC-CO point-to-point correlation at 2" 

In order to compare point-to-point RG, GO, and mid-IR bright- 
nesses in the observed regions, we overlaid on the images regu- 
lar grids of rectangular boxes, whose dimensions coincide with 
the beam of observations (corresponding to a spatial scale of 
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X 


y 
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r 


slope 
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r 


Ico 


Irc 


0.43 


0.17 


0.56 


0.77 


0.21 


0.84 


Irc 




1.47 


0.22 


0.69 


1.11 


0.25 


0.74 


Irc 




1.19 


0.16 


0.85 


1.02 


0.13 


0.95 


Irc 


Is/zm 


0.99 


0.14 


0.83 


0.87 


0.14 


0.84 


Ico 


l4.5/im 


0.52 


0.3 


0.41 


0.87 


0.34 


0.53 


Ico 


Is.S/iiii 


0.49 


0.21 


0.40 


0.78 


0.25 


0.77 


Ico 


Ig/zm 


0.39 


0.17 


0.52 


0.69 


0.23 


0.67 



Table 3. Coefficients of IR-RC-CO 2"correlations. Observed 
data have been fitted with the function y(x) = c ■ x''°'"'. Slope 
is the result of the fit, rms is expressed in dex, and r is the 
Pearson's correlation coefficient. 



100 pc), and we averaged all pixel values within each box. We 
performed the power law fitting described above, determining 
for each correlation the corresponding slope and normalization 
parameter. The results of this fit are reported as labels on the 
plots. 

Figures|7]and[8]show the overlay of RC brightness contours 
on CO integrated intensity images and the plot of 1 1.4 versus Ico 
for both regions. The overlay of the emissions shows that CO 
and RC emissions distribute very similarly in both the north 
and south bar ends, despite the presence of minor local dif- 
ferences. In particular, both RC and CO emissions follow the 
same filamentary distribution in both regions. The plot of the 
point-to-point correlations show that in the north bar end the 
RC emission spans only one order of magnitude, whereas the 
CO emission varies over two orders of magnitude. In the south 
bar end instead, both RC and CO emissions vary by a compa- 
rable factor This results in a flatter slope in the northern region 
(0.43) than in the southern one (0.77). The Pearson's correla- 
tion coefficient, r, is also different in the two regions (r=0.56 
in the north bar end and r=0.84 in the south bar end), indi- 
cating that the RC is less correlated to the CO in the northern 
region than in the southern one. The scatter of both correlations 
(0. 17 dex and 0.21 dex in the northern and southern region, re- 
spectively) is similar to that obse rved in the 6" correlations 
presented in lPaladino et al. I (l2006h . 



We performed the same point-to-point analysis using the 
4.5, 5.8, and 8 jum "dust maps" (see Sect. 12.31 ), convolved to 
the same beam of RC and CO images. The results for north 
and south bar ends are reported in Figs. |9] and [TOl respectively. 
In the two regions, it appears that both the RC and mid-lR 
emissions vary by about one order of magnitude, whereas the 
CO varies over almost three orders of magnitude, producing 
flat CO-mid IR correlations. The RC is well correlated with 
dust emissions of both 5.8 and 8 yum: the resulting correlations 
are nearly linear and scatters vary from 0.16 to 0.13 dex. The 
RC-4.5 jum correlation has a higher scatter, in particular in the 
southern region. We report, in Tabled the fitted slopes, rms of 
the 2" correlations, and the Pearson's coefficients. 



3.2. The comparison with Ha images and the 
determination of thermal radio flux 

The two regions of NGC 3627 show very different stellar and 
dust distributions. Some wor ks report on o bser vations of HII 



regio ns in this galaxy (see iHodgd 1 19741 and IChemin et al 



20031) . In the north bar end, crowded with stellar objects, dust 
and HII regions have been the subject of various works focused 
on s upernovae SN1989B and SN1973R ho sted in NGC 3627 
re.g.. lVan Dvk et al.lll999[rWells et al.lll994 . 

Figure [TT] shows the comparison between the RC and CO 
emissions (contours) an d the po s itions of brightest Hll regions 
(crosses), as reported in lHoded (Il974l) . Some distinct sources 
of RC emission lie over the positions of HII regions. Without a 
radio spectral index to verify the thermal nature of the radiation 
from these sources, we can only tentatively identify the radio 
emission as originating in the HII regions. Indeed, superposing 
the north bar end RC and Ha images, we found that the spatial 
distribution of the brightest Ha knots is well matched by the 
RC peaks of emission. Instead, the CO emission is located on 
the edges of Hll regions, tracing the dust arc-like distribution. 

In the south bar end region the RC and CO emissions show 
an extended distribution similar to that of dust, without re- 
solved peaks of emission. In the Ha i mage, there ar e some 
distinct bright knots also identified by iHodgd (119741) as HII 
regions, but they do not have a point-like RC-counterpart. 

We used the H^^ image to estimate the RC thermal com- 
ponent to evaluate its eff'ect on the correlations. The typical 
thermal fraction of radio emission at 1 GHz, measured from 
integrated spectra of spiral galax ies is of the o rder of 0.10. 
In particular, for NGC 3627, Nik las etal.lll997l determined a 
thermal fraction < 0.09. We estimated the thermal radio flux 
from the emission i n the Ha recombination line (see, e.g., 
Cjplan & Deharven j [l986l and iNiklas et al.iri997h . using flie 
iNiklas et all (11997!) Eq. 2 in the following way: 



= 2.238 -lO^ ■ 

VmJy / 



I Sua \ /Te\0-42 

Vergs^'cm^^/ \K/ 



0.04995 /Te\ 



(1) 



where Sth is the thermal radio flux density at frequency v=1.4 
GHz, assuming a thermal electron temperature Tg of lO'* K and 
measuring the Ha flux: Sho- The mean value of the thermal 
flux, measured in the same grid of point-to-point analysis in 
the northern region, is 0.02 mJy. This value rises to 0.08 mJy in 
correspondence with the brightest Hll regions (??, indicated as 
bold-face crosses in Fig.fTTT). In the northern region, we found 
that the thermal fraction ranges from 0.07 to 0. 11 . Even if our 
measure underestimates the thermal emission, owing to absorp- 
tion of the Ha emission by dust, it is still likely sufficient for 
the purpose of roughly estimating the thermal flux density and 
to ensure that the thermal contribution to the radio emission on 
the spatial scales analyzed is of order of 10%, not significantly 



differ ent from that measured on integrated spectra ( INiklas et al 
1997h . 
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S/jm-RC-CO 



name 


4> 


e 


rms 




e 


rms 


6" correlations 


NGC 3351 


0.94 


0.61 


0.16 


0.77 


0.69 


0.13 


NGC 3521 


0.77 


0.35 


0.14 


0.57 


0.39 


0.11 


NGC 3627 


0.95 


0.60 


0.23 


0.63 


0.75 


0.19 


NGC 4826 


0.92 


0.74 


0.15 


0.76 


0.82 


0.14 


NGC 5 194 


0.97 


0.52 


0.23 


0.81 


0.60 


0.20 


NGC 7331 


0.87 


0.36 


0.10 


0.73 


0.38 


0.09 


All 


0.92 


0.60 


0.25 


0.65 


0.72 


0.24 


2" correlations 


NGC 3627 N 
NGC 3627 S 








0.32 
0.58 


0.32 
0.62 


0.19 
0.23 



Table 4. 3D fit coefficients. Observed data have been fitted with 
a line in three-dimensional space. and (p are the direction an- 
gles of the line, rms is the scatter of the data. 



3.3. The 3D view 

The RC-CO-IR correlations are usually studied separately as 
projection in different planes. What we observe is indeed a cor- 
relation between all three emissions. This three-dimensional 
aspect cannot show up in the usual 2-dimensional description. 
Hence, we present a new perspective, a three-dimensional rep- 
resentation, with which to visualize and study the RC-CO- 
IR correlation. We performed a fit to the data in a three- 
dimensional space, placing log(Ico)^ log(IiR), and log(lRc) in 
the X, y and z axis, respectively. We computed the line fitting of 
the observed points, by minimizing the distance of these points 
from the line: 

{X — xq + p ■ cos6cos(p 
y = yo + P ■ cos9 sin(p 
z-zo+p - sinO 



where {xQ,yQ,ZQ) is a generic point through which the line 
passes, 9 and (f) are its direction angles, and p is a parame- 
ter. As the fit procedures are completely different, the 3D fit 
coefficients cannot be easily related to the power law indexes 
found for the two-dimensional correlations. However, our pur- 
pose is to introduce a promising technique to visualize the "uni- 
fied" correlation between the three emissions. The extension of 
this representation to a larger sample of galaxies could possi- 
bly give insight into the physical explanation of the correlation 
itself. 

We used the six galaxies studied in I Paladino et al. (2006) 
as a reference sample for the three-dimensional correlation. 
Figure [12] shows the RC-CO-24 jjm, with the plotted line rep- 
resenting the fit to all the data. We conducted the same analysis 
on RC-CO-8 fim correlation, but for the sake of conciseness, 
this figure is not shown. Values for the angles and (p, and the 
rms, resulting from the fit, for each galaxy are reported in Table 



the scatter among different galaxies, whereas the scatter within 
each galaxy is even smaller (see Table|4]l. 

Figure [T3] shows the comparison between the three- 
dimensional 2" RC-CO-8 im\ correlation in the two observed 
regions and the 6" correlation on the entire galaxy. Red boxes 
and blue triangles represent measures obtained in the north and 
south bar ends, respectively. The fitted lines for each set of data 
are also plotted. The bold line is the result of the fit to the data 
of the entire galaxy at 6". The line fitting data of the south 
bar end region presents direction angles similar to those of the 
entire galaxy, and resembling to those of all galaxies. On the 
contrary, the line fitting the northern region data has lower di- 
rection angles, indicating a flatter correlation. 

4. Discussion 

Studying the RC, CO, and 8 jim emissions from two regions 
of spiral galaxy NGC 3627, we found that, despite some mi- 
nor local variations, the spatial distributions of RC, CO, and 8 
/zm emissions match very well. The phenomeno logical model 



(see | Bicay & Heloui 11990. and, more recently. Murphy et al 



In these correlations, we deal with surface brightnesses, and 
thus they are unaffected by distance biases. The scatter of the 
global 3-D correlation is of 0.24 dex, i.e., a factor 1.7. This is 



! 2006b|) proposed that the radio images should be a smoothed 
version of the IR maps due to the diffusion of cosmic-ray elec- 
trons from the star-forming regions. The results obtained in our 
study may suggest a more complex picture. 

The propagation of cosmic rays along magnetic field lines 
is, however, still matter of debate. It is generally recognized 
that the cosmic rays propagating along the galactic magnetic 
field are continuously scattered and accelerated (see iMelrose , 
1968b . Since the late 1970's, the proposed theory of diffusive 
shock acceleration has been developed and widely accepted. 
This theory proposes the statistical Fermi's acceleration mech- 
anism as the responsible for the particles acceleration (IFermi 
1949b and charges the slowing of such particles with reflec- 
tions at strong sh ock waves. Without going into detail (see 
Strong et al.ll2007[ for a recent review on the cosmic ray prop- 
agation and interactions in the Galaxy) we can say that the 
bulk propagation of cosmic rays is slowed down to the velocity 
of magneto-hydrodynamic (MHD) waves, which is the Alfven 
speed, Va = ^JB^/87^p = 1.6 ■ 10^ n"'^^ B^g, in cm s"', where 
B is the magnetic field strength, and p and n are the density 
and number density of the gas, re spectively. For a typical disk 
magnetic field strength B~ 5 fiG (ICondonll 1992b and a density 
of the warm phase of the galactic disk ISM (which may be re- 
sponsible fo£_turbulent diffusion of cosmic rays) n=0.3 cm"-^ 
(iBiskamp l2003b . the Alfven velocity is ~15 km s '. The syn- 
chrotron lifetime t, at Vc=l GHz, is of about 10^ yrs. In these 
conditions, the diffusion scale-length of the radiating electrons, 
Ldiif ~ Va ■ T, is of 1.5 kpc. 

The scale we have examined (100 pc) is much lower than 
the electron diffusion scale. Nevertheless, in the two observed 
regions the RC emission does not appear to be more diffuse 
than that of CO and 8 /vm. We compare the relative shapes of 
the emission profiles taking a one-dimensional slice of the RC, 
CO, and IR images. Due to the clumpy distribution of the emis- 
sions at this resolution, the profiles obtained along some paths 
contain multiple emission maxima corresponding to resolved 
features. In Fig. [141 the profiles obtained along two slices in 
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SOUTH B,\H END 



NORTH B.\H END 





on ph oto-dissociation regions see, e.g., iHoUenbach & Tielens 



Fig. 14. Emission profiles of RC, CO and IR. Left panel: 
Profiles along the horizontal direction in the south bar end re- 
gion. Right panel: Profiles along a slice in the thin arm. See 
Fig.[TT]for the slice position. 



1999h . This option can explain the lack of CO emission in 



correspondence with bright Hll regions with RC counterpart, 
nevertheless it still need to be quantitatively verified. Certainly, 
more complex mechanisms are taking place in the observed re- 
gions, especially since the ends of the bar might suffer from 
turbulent motions that influence the velocity field. 

Since the RC brightness critically depends on the galaxy's 
magnetic field B and cosmic -ray electron distribution No, and 
the separation between these two contributions is not yet clear, 
the interpretation of the results obtained is difficult. Spatially 
resolved spectral index studies at low frequencies, where the 
thermal emission component is seemingly negligible, can help 
to clarify the situation. 



the observed regions (the positions of slices are shown as black 
lines in Fig. fTTT i are reported as an example. The IR, RC, and 
CO profiles have been plotted with normalized peak amplitude 
to facilitate a comparison of their shapes and their spatial co- 
incidence. The slice in the southern region is taken along the 
horizontal direction. The RC and CO widths are comparable 
and not greater, as possibly expected, than the analogous width 
of the dust emission at 8 //m. For the northern region we re- 
port the profiles taken along a slice in the thin arm structure, 
where all the three emission are detected. In this case the most 
concentrated emission appears to be the CO one, while the RC 
emission profile does not appear to be a smoothed version of 
the IR one. In both regions, there is evidence that RC emission 
is not more diffuse than the CO and 8 //m emission. This is a 
surprising result, as the observed regions are much smaller than 
the expected electron diffusion scale. 

It is important to stress that we found that the RC emission 
is well correlated to the emissions at yu m wavelengths in the 
northern and southern region. The yum emission traces a dif- 
ferent dust component than that traced by the FIR emission. 
In practic, while the latter represents cold dust, as mentioned 
before, the 4.5 jjm emission contains component of hot dust, 
whereas in 5.8 and 8;um images the emission from PAH be- 
comes significant. The CO emission varies over three order of 
magnitude while dust emission spans only one order of magni- 
tude. In our opinion, this result could be explained considering 
that the CO emission traces molecular gas where star forma- 
tion may just begin, whereas the bulk of the 8 p.m. emission, for 
example, comes from PAH molecules and these particles are in 
turn concentrated near regions of current star formation. 

As we noted in previous sections, the north bar end region 
presents flat correlations of CO emission with both RC and IR. 
The flat RC-CO slope can be produced either by an higher RC 
brightness or a lower CO integrated intensity with respect to 
the average RC-CO relation. Assuming that the excess of RC 
brightness is due to the RC thermal component, the thermal 
fraction would be of 0.55. This value is higher than what we 
measure from the Hq. emission (see Sect. 13.2b . and thus this hy- 
pothesis is unlikely. Otherwise, if the CO is lower than average, 
a possible explanation can be the photo-destruction of molecu- 
lar clouds near HII regions by ultraviolet photons (for a review 



5. Summary and concluding remarks 

We studied the correlations between the 1.4 GHz radio con- 
tinuum, the CO, and the IR emission at short-wavelengths in 
two regions of the spkal galaxy NGC 3627. We carried out new 
VLA and PdBI observations down to an angular resolution of 
order of 2", which corresponds to a resolution of ~ 100 pc at 
the distance of 1 1 Mpc. We used the released IRAC and Ha 
SINGS images of comparable resolution. Even if the scale we 
examined is much lower than the electron diffusion scale we 
do not find evidence of a diffusion of the RC emission with re- 
spect to the scales of CO and 8 yum, as would be expected by 
previous phenomenological models. 

The observed regions are located at the northern and south- 
ern ends of the bar of NGC 3627 where two broad peaks of 
emission are shown both in RC and CO images. Our new high 
resolution observations permit a distinction between compact 
and extended emission that would not be possible with previ- 
ous data. In particular, in the northern region we can resolve 
some RC peaks have a spatial coincidence with HII knots. 

We point out that our results refer to one single object, and, 
in particular, to two distinct regions in this object, and we do 
not attempt to provide a definite answer to the question. Our 
main result is that in these two regions, even at scales of 100 pc, 
where we would expect a break down of the correlations due to 
the cosmic rays diffusion, the RC emission is still correlated to 
that of CO. Furthermore we found that the RC emission is also 
correlated to the 4.5, 5.8 and 8yL(m dust emission. Whatever 
these emissions trace, their correlation with RC emission on 
such a small scales is surprising and requires explaination. The 
study of a significant sample of galaxies would answer some of 
the questions left open. 

Interestingly, in the northern region the RC-CO and IR- 
CO correlations have flatter slopes with respect to the aver- 
age correlations, which might be due, among other options, to 
the photo-dissociation of molecular clouds near HII regions. 
Despite this minor deviation, the spatial distributions of RC fol- 
lows well those of CO and IR, tracing thin arc-like structures. 
The widths of several structures, measured along slices across 
them, are similar in the RC, CO, and 8 ;um maps. This indicates 
that the scale of RC emission is not larger than those of 8 jim 
and CO. 
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Furthermore, to give a complete description of the RC-CO, 
CO-mid IR, and RC-mid IR correlations, we have proposed 
a more concise and effective three-dimensional representation 
of the correlation, producing a three-dimensional fit of the ob- 
served data. We found that RC, CO, and mid-lR observed data 
are fitted by a line in three-dimensional space with a scatter of 
about 0.25 dex. In our opinion, this is a promising technique to 
visualize the "unified" correlation between the three emissions. 

This study is focused on two individual regions of the 
galaxy NGC 3627, and the analysis of more objects has cer- 
tainly to be done. Nevertheless, the observed correlation be- 
tween RC, CO, and mid-lR emissions down to spatial scale 
of 100 pc may suggest that the coupling between electrons 
diffusion, losses, and injection is far more complex than ex- 
pected. Future works on spectral index at high resolution will 
help to separate the contributions to the RC brightness due to 
the strength of magnetic field from that due to the density of 
relativistic electrons. 

Appendix A. C0(2-1) observations 

Figures [15] and [16] show the '^CO(2-l) velocity-channel maps 
and the integrated intensity observed with the PdBI in the north 
and south bar end regions, respectively. 

The '2CO(2-l) and '2CO(1-0) emissions are detected in al- 
most the same velocity range. The 230 GHz channels have been 
compressed to have the velocity resolution of 3.26 Km s"'. 

Due to the particularly severe missing flux problem in 
these images, they are not reliable for quantitative studies. 
Nevertheless, we note that the northern '^CO(2-l) emission 
has a peak coincident with the '^CO(l-O) peak, and that in the 
southern region three different peaks of emission are resolved 
(Figure [T6]l. 
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Fig. 7. Left: Overlay of RC brightness (contours) on CO integrated intensity (grey-scale image) for the north bar region observed 
at a resolution of 2 '.'74x1 '.'76. Contours level start from 0.2 mJy beam'' (5<t) and scale by a factor of V2. Note that the spatial 
coincidence between the two emissions is maintained well below the kpc scale. Right: I1.4 versus Ico- Brightnesses are measured 
by averaging all pixel values within beam-sized boxes recovering the entke region of emission. Only points above 3 cr level have 
been plotted. 
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Fig. 8. Left: Overlay of RC brightness (contours) on CO integrated intensity (grey-scale image) for the south bar region observed 
at resolution of 2'.'74xl '.'76. Contours level start from 0.15 mJy beam"' (5cr) and scale by a factor of V2. Note that the spatial 
coincidence between the two emissions is maintained well below the kpc scale. Right: Ii 4 versus Ico- Brightnesses are measured 
by averaging all pixel values within beam-sized boxes recovering entire region of emission. Only points above 3 cr level have 
been plotted. 
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Fig. 11. Comparison between radio continuum (upper panels) and CO (lower panels) maps with HII regions positions. RC and 
CO contours start from 0.2 mJy beam ' and 2 Jy beam ' km s respectively, and scale by a factor of 2, to enhance point 
sources. Red crosses represent positions of br ightest HII reg ions in the north (left panels) end, south (right panel) bar end. For a 
complete list of HII regions in NGC 3627, see iHodgd (11974 . Black lines in upper panels show the positions of two slices for the 
comparison of emission profiles (see Sect. 4). 
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Fig. 12. RC, CO, and 24 fim 3D correlation. The x, y, and z axis are the log(CO), log(IR) (24 //m), and log(RC), respectively. 
Data at 6"resolution of six galaxies reported in Table |4] are shown and the line plotted is the result of the fit for all data. The 
scatter of this fit is 0.25 dex. 
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Fig. 13. RC, CO and 8 fim 3D correlation at 2". The x, y, and z axis are log(CO), log(IR) (8 /urn), and log(RC), respectively. 
Data at ~ 2"resolution of northern and southern region of NGC 3627 are represented by red crosses and blue triangles, respec- 
tively. Red solid and blue dotted lines represent the fit of these two set of data. The bolded black line corresponds to the fit of 
6"RC, CO, and 8 urn data of the overall galaxy NGC 3627. 
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Fig. 15. North bar end of NGC 3627. Upper panel: '^CO(2-l) velocity-channel maps observed with the PdBI. The spatial 
resolution is of 1 '.'59xO'.'88 at PA=22°. We map an area of view of 15". The phase tracking center is indicated by a cross at 
Q'y2000=ll''20'"13.5' and 5y200o=13°00'17'.'7. Velocity-channels are displayed from v=-53.76 km s"' to v=-148.26 km s"' in 
steps of 3.26 km s Velocities are referred to theLSR scale and the zero velocity corresponds to vo=712.6km s Contour levels 
start from -20 mJy beam"' with step of 20 mJy beam"'. The rms noise is 3.4 mJy beam"' and only regions whose brightness is 
larger than 6-cr are shown. 

Lower panel: '^CO(2-l) integrated intensity. Contours level start from 0. 18 Jy km s"' beam"' and scale by a factor of V2. 
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Fig. 16. South bar end of NGC 3627. Upper panel: ^'CO(2-l) velocity-channel maps observed with the PdBl. The spatial 
resolution is of 1 '.'63xO'.'87 at PA=31°. We show a field of view of 15". The phase tracking center is indicated by a cross at 
Q'/2000=l 1''20'" 16.6' and 5y2()00=12°58'44'.'5. Velocity-channels are displayed from v=216.72km s^' to v=92.89 km s^^ in steps 
of 3.26 km s Velocities are referred to the LSR and the zero velocity corresponds to vo=712.6 km s ^ Contour levels start 
from -20 mJy beam"' with step of 20 mJy beam"'. The rms noise is of 4 mJy beam"' and only regions whose brightness is 
larger than 5-cr are shown. 

Lower panel: South bar end of NGC 3627:'^CO(2-l) integrated intensity. Contours level start from 0.24 Jy km s"' beam"' and 
scale by a factor of V2. 



